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ABSTRACT 

We report on the status of quark and gluon propagators in quenched, gaugefixed lattice QCD. In 
Landau gauge we find that the effective quark mass in the chiral limit is M q ~ 350(40)MeV. 
Quark and gluon propagators, the slope of the quark dispersion relation, and effective masses all 
appear to depend on gauge. A link-chain picture of lattice gaugefixing in the color N — > OO and 
strong coupling limit, where the system becomes almost solvable, supports the gauge variance of 
these numerical results. 

1. Motivation and Overview 

Lattice gaugefixing has drawn considerable attention in recent years. Gluon [1], 
quark and photon propagators and effective masses have been evaluated on the lattice 
in Landau, Coulomb and related "A-gauges, " A<9 A + V • A = 0. Wavefunctions of 
valence quarks inside mesons have been probed with gauge- variant extended sources. 
Furthermore, lattice gaugefixing is routinely performed nowadays in order to use 
extended gauge- variant sources as interpolating fields in the calculation of gauge- 
m variant quantities. While beyond the scope of this talk (see Refs. [2]- [4] for some 
references), the numerical work has spawned analytical and computational studies 
of longitudinal and topological gaugefixing ambiguities and their effects on gluon, 
quark and photon correlation functions and operator product expansion coefficients 
determined from gauge covariant matching conditions. 

For perspective let us focus momentarily on the solvable Schwinger model, where 
quarks(electrons) are confined but the photon is physical and has a gauge-invariant 
nonzero mass from the £7a(1) anomaly. Coulomb gauge — the unitary gauge for the 
Schwinger model — has no unphysical modes and the Coulomb gauge quark propagator 
is the physical amplitude for quark propagation. Since in f + f dimensional empty 
space a quark's electric field E is constant and cannot die away, E produces a constant 
flux of quark-antiquark pairs moving away from the original quark at the speed of 
light. Consequently the quark propagator — the amplitude to have just one quark 
at x ^ starting with an x = quark — vanishes and the effective Coulomb gauge 
quark mass diverges. On the other hand, the covariant gauge quark propagator is not 
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the physical amplitude for quark propagation. In fact the Landau gauge quark mass 
vanishes. This discrepancy with Coulomb gauge is due to unphysical modes (or, 
alternatively, Gupta-Bleuler physical-state conditions on the Hilbert space) which 
ruin the identification of the covariant gauge quark propagator with the physical 
propagation amplitude. Thus, the effective quark mass can be shown to vary with 
covariant gauge parameter [4]. The photon mass, being physical, is gauge m variant. 

The point is, in association with confinement, the effective quark mass is ambigu- 
ous because it varies with gauge — at least in the Schwinger model. Like Schwinger 
model quarks, QCD quarks and gluons are confined and unphysical and nothing guar- 
antees their masses be gauge-invariant. If what is not forbidden occurs (as usually 
happens), quark and gluon masses would be gauge- variant. 

2. Numerical Lattice QCD Results 

Numerical quark and gluon propagators at large times — as large as our lattices 
permit — look approximately, but not exactly, like free particle propagators. 2 In A- 
gauges M q and M g are evaluated by matching gaugefixed p = propagators to 
free particle propagators at large Euclidean time t. In particular, the zero mo- 
mentum A-gauge Af g ( A ) is obtained by comparing the numerical quark propagator to 

Z q ( 1+ 2 )e~ M v ) *. The "(A)" superscript allows for A dependence. Since the scalar "I" 
part of the numerical quark propagator equals the vector "7 " part at large enough 
times, we can focus on the scalar part without loss of generality. Correlated x 2 -fits of 
the large time scalar part of lattice quark propagators to Z q e~ M i * give unacceptably 
large % 2 s. Quark propagators fit well to trial functions like At B e~ mBt . 

Rather than use an unmotivated but better-fitting trial function, in the rest of this 
talk we will ignore the large x 2 's and assume free particle form Z q e~ M i *. Then the 
following features are observed: (i)The Dirac scalar component of the p = numerical 
quark propagator equals the vector component at large times. (ii)As depicted in 
Figure 1, quark and gluon effective masses vary with gauge parameter A. They 
decrease with increasing A such that most of the change occurs when < A < 2. 
As A — > 00, the masses stabilize to a nonzero value. (iii)Motivated by the chiral 
perturbation theory relation Af 2 oc m q , a fit to 

Mf) = 6 (A) M 2 + M C W (I) 

in the chiral regime yields in Landau gauge at /3 = 6.0 that &W ~ 2.7(.3) x I0" 4 /MeV 
and MP ~ 350(40) MeV. In general, M q is about one-half the p mass at the same 
hopping constant K, the lattice bare mass parameter, independently of K. (iv)M g is 
about 70% larger than . (v)The slope of the dependence of quark energy square, 
-E 2 , on spatial momentum square, p 2 , varies with A so that not all gauges can be 
consistent with Lorentz invariance. 

2 The numerical calculation is on (3 = 5.7 and 6.0 quenched Wilson lattices. For noncogniscenti, 
[3 is a measure of physical lattice spacing; the larger [3 is, the smaller the lattice spacing and the 
closer one is to the continuum limit. (3 = 6.0 correspond to a lattice spacing of ~ l.TGeV -1 . 
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FIG. 1. Comparison of to Strong Coupling A-Dependence 



3. Lattice Gaugefixing and f3 = 0, N — >■ oo Model 

The A dependences of the effective quark mass and dispersion relation are quali- 
tatively reproduced in a f3 = 0, color N — > oo analytical solution of A-gauge lattice 
QCD. Figure I compares the chiral quark mass at strong coupling to the nu- 

merical effective quark masses rescaled so that they all agree at A = I . The strong 
coupling calculation was done in the following way. A-gauge vacuum expectation 
value (0) in lattice QCD is defined as 

(0) EE [[[©]/]J u , [©], EE ^/[^e-^0, Zg EE /[^'K^ (2) 

where labels 6 = /, v or u refer respectively to quark, gauge transformation, and 
lattice gauge "link" fields. [0]# is the average of the operator (or combination of 
operators) over the fields 6 = /, v, or u. f stands for quark fields ip; v } for fields V 
which transform a background gauge field into A-gauge; and it, for the lattice gauge 
fields U. S-f and S u are the usual lattice Wilson fermion and plaquette actions. The 
lattice gaugefixing action S v is 

SV = -J2'Y Retr(K^yJ +A ), J, ee ^ (3) 

where A ee A and A; ee 1. Quarks are quenched by choice; consistency with the 
Fadeev-Popov method requires {V x } to be quenched. S v is minimized site-by-site 
with respect to {V x } in numerical simulations to achieve A-gauge. Similarly, the 
£ — > limit corresponds to A-gauge. 



Since this system is not solvable, we take the A-gauge, N — > oo limit [3]. Identi- 
fying £ = l/(2iV), taking N — > oo and integrating out gaugefixing fields {V x } yields 

(V y (yA w )vj) a] = 8 aj £ A w » y ( J ) tr ^- A w » y = -^ tT ( v yyK v h\fJ=o- (4) 

a and j are color indices ranging from 1 to N and y A w is a link chain from w to 
y. Eq. (4) says the gaugefixed expectation value of y A w is the weighted sum over all 
Wilson loops made by joining to y A w a self avoiding link segment w A y . (4) identifies 
weights A j with the f3 = gaugefixed expectation value of w A y . All of the A 
dependence is in these weights. 

Hopping expanding the quark fields, applying (4), and resumming the hopping 
expansion leads to an expression for the quark propagator. In D = 3 + 1 dimensions 
in the color N — > oo and f3 = limit, and in an orthogonal trace approximation 
described in Ref. [3], the quark propagator dispersion relation is 



£(A)2 = M (A) 2 + E ^2 5 g i{X ) = !_l L (5) 



Since slope g 8 (A ^ 1) ^ 1, the dispersion relation is not free-particle form except in 
Landau gauge. At p = ~g i dependence drops out, but remains A-dependent. 

The quark propagator in this limit (with Wilson parameter r = 1) has the free particle 
exponential form. 

While we do not explain Wilson parameter r here, let us note that in strong 
coupling the r = point (where backtracking and hence spontaneous chiral symmetry 
breaking is allowed) is more similar than r = 1 to continuum QCD. At r = 
varies with A, diverges in Coulomb gauge, and is analytic and linear in m q near 
m q — > 0. Following (1), the Taylor expansion = + B^m q + 0{m\) yields 

M (A) _ { 67 - 25X2 5 (A) _/ 2Jm ^^ A<§; 

1V1 c — 0A0 \ ) -25 + 200A+672A 2 \ > 1. D ~ \ 4753-38024A+399296A 2 \ \ 1 l D J 
Z ^ ZA I 4(4A-1) — 2 ' [ 117128A(4A-1) A - 2" 
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